Exact Results for the Reactivity of a Single-File System 
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We derive analytical expressions for the reactivity of a Single-File System with fast diffusion and 
adsorption and desorption at one end. If the conversion reaction is fast, then the reactivity depends 
only very weakly on the system size, and the conversion is about 100%. If the reaction is slow, then 
the reactivity becomes proportional to the system size, the loading, and the reaction rate constant. 
If the system size increases the reactivity goes to the geometric mean of the reaction rate constant 
and the rate of adsorption and desorption. For large systems the number of nonconverted particles 
decreases exponentially with distance from the adsorption/desorption end. 

PACS numbers: 02.50.Ga, 68.65.-k, 82.75.Qt 



Zeolitic structures like mordenite, ZSM-22, and 
AIPO4-05 are industrially immensely important. They 
are examples of single-file systems (SFS).0 They have 
one-dimensional pores through which molecules can dif- 
fuse, but with cross sections that are too small to allow 
passing. As for many other ID systems, this leads to very 
interesting kinetic effects. |Q Research on SFS zeolites has 
focussed on the diffusion, because the mean-square dis- 
placement in a SFS is proportional to the square root 
of time, and not proportional to time.[^, ^, ^ Molecular 
Dynamics, I, |[ § Dynamic Monte Carlo (DMC),|, |, 
[l0| and reaction-diffusion equations have mainly been 
used. Few studies have included reactions, jl^ [l^, |l3| Wc 
have argued that this is unfortunate, because if diffusion 
is fast the temporal dependence of the mean-square dis- 
placement is not very relevant, but there are still effects 
from the non-passing of the particles in the pores. 

In a previous paper we have investigated steady-state 
properties of a SFS for different assumptions of the re- 
active site distribution using DMC as well as analytical 
techniques. jlOt In the current paper we derive exact re- 
sults for the case that the entire pore is reactive. These 
results give detailed insight in the relationships of the 
various system parameters and the reactivity. 

Our model consists of S sites forming a one- 
dimensional finite chain. Each site is numbered consec- 
utively from 1 for the site on one end, to S for the site 
on the other end. Each site is vacant or is occupied by a 
particle. We have two types of particles: A and B. An A 
can be converted into a B on any site. Adsorption and 
desorption of particles can only occur on site I. Only A's 
are adsorbed, but both A's and B's desorb. Both types 
of particle diffuse by making random hops to neighboring 
sites if vacant. 

The evolution of the system is described by a master 
equation |l^ 



dPg 
dt 



J2 \Wo.pPp - WpM , 
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(a particular distribution of particle over the sites), the 
P's are the probabilities of the configurations, t is real 
time, and the W's are constants that give the rates with 
which reactions change the occupations of the sites. Wap 
corresponds to the reaction that changes /3 into a. The 
rate constants in our model are Wads for adsorption of an 
A onto site I if vacant, Wdos for desorption of a particle 
from site 1, Wrx for the conversion of an A into a B on 
any site, and Wdis for a hop of a particle to a vacant 
neighboring site. 

If the diffusion is infinitely fast, we can derive a simpler 
master equation for the number of particles in the system. 
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dt 
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+ Wdes^—PN + l-Wdcs^PN 
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where a and (3 refer to the configuration of the adlayer 



where Pn is the probability that there are N particles 
in the system. This is a master equation of a one-step 
Markov process. |l5| 

We are interested in the probability distribution 
fMK{t) that if at time t = the number of particles 
is K, this number becomes M for the first time at time t 
with M < K. If there are N particles in the system and 
at time i = a particle adsorbs, then fN,N+i is the prob- 
ability distribution for the time that this particle desorbs. 
This is based on the property of a SFS that particles can 
not pass each other. So /jv,jv+i is the probability distri- 
bution for the residence time of a particle that adsorbs 
in a system with N particles. 

Let Pnk be the solution of Eq. (||) with TV the number 
of particles and the initial condition Pata'IO) = Snk- Let 
Qnk also be a solution with Qnk{0) = ^nk^ but now for 
the master equation with an adsorbing boundary at M: 
i.e., we remove the term in Eq. (j^) that corresponds to 
an adsorption then there are M particles in the system. 
With N > M we have 

PNK(t) = QNKit) + I dt'PNAlit - t')fMK{t'). (3) 
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This equation is called the renewal equation. If we 
take N = M in the renewal equation, then we have 
Qhixit) = by definition. So we get the integral equa- 
tion 



PMK{t)^ / dt'PMM{t-t')fMK{t') 



(4) 



for fMK- 

A particle that adsorbs, stays in the system for a period 
t, and then desorbs, has a probability exp(— Wrxi) that 
it is not converted during the time it was in the system. 
Because the probability distribution for the time that it 
stayed in the system is given by fN,N+i{t), the proba- 
bility that the particle desorbs without being converted 
is 



dt e^p{-Wr^t)fN,N+i{t)- 



(5) 



This is equal to fN,N+i{Wrx) with fN,N+i the Laplace 
transform of fN,N+i- The reason for using this Laplace 
transform is that is that it is related to the Laplace trans- 
form of the solution of the master equation through the 
renewal equation. Laplace transforming Eq. @ yields 



Pmk{s) = PMM{s)fMK{s)- 



(6) 



If we write the master equation (g) in matrix-vector no- 
tation as P = WP, and Laplace transform it we get 



Y,{s5nm - ^nm)Pm{s) = Pjv(O) 

M 

With the initial condition for Pmk this yields 
Pmk{s) - [(s- W)-ii 



so that 



fMK{s) 



[{s - W) 



MK , 



MK 
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(7) 
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We define the reactivity Pprod sls the number of par- 
ticles that is being converted per unit time. To get 
Pprod we have to multiply the probability that a par- 
ticle is converted by the rate of adsorption and make 
a weighted average over the number of particles in the 
system. The probability that a particle is converted 
equals 1 — /Ar,jv+i(W^rx), and the rate of adsorption equals 
W^ads(l — N/S). At steady state the number of particles 
in the system is given by ||l0|] 



P; 



{ss) 



N 





S 


rvF-adsi 


N 











(10) 



Combining this results in 



s-i 

N=0 



N 



s 



N,N+l{Wrx 



(11) 
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FIG. 1: Reactivity as a function of rate of adsorption and 
desorption. The reactivity and the rate of adsorption and 
desorption are scaled with the rate constant of the conversion. 



Figure]^ shows this reactivity for two finite system sizes 
and three loadings. The loading 9 is defined as the prob- 
ability that a site is occupied at steady state. It's equal 
tofBl 



ads 



W^ds + Wdc 



(12) 



The reactivity is compared to the (average) rate of ad- 
sorption and desorption A, which at steady state is given 
by 



A = I^ads(l -0)^ Wdc 



(13) 



If the rates of adsorption and desorption are small com- 
pared to the reaction rate constant Wix, then Ppiod — * 
This means that almost every particle that enters the 
system is converted. If the rates of adsorption and des- 
orption are much larger, then 



B 



prod 



(14) 



for A/Wrx oo and S not too large. Figure |^ shows 
the 9 dependence for small S and the S dependence for 
small 9. It does not show this S and 9 dependence for 
S and 9 both large. A needs to be much larger than the 
values in the figure for that. For S ~ 100 we see that 
the reactivity may be proportional to 9 for 9 < 0.5, but 
not around 9 — 0.9. For the latter value the reactivity is 
even smaller than for 6 = 0.5. 

In practice only Wads can usually be varied indepen- 
dently, either through the pressure or the concentration 
of A's outside the system. The limit W^ads ^ oo corre- 
sponds to 6* = 1 and A = Wdos- So we see that there 
generally is an upper bound on A. 
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Comparing S = 5 and S = 100 in Fig. |i| for some finite 
A shows that iJprod varies less with 6 for larger S. We 
can show that in the limit S ^ oo the reactivity depends 
only on A and W^^. If the system becomes infinitely large, 
then the fluctuations become small with respect to the 
number of sites. If the system is then at steady state we 
can write the master equation as 



N 



dt 



XiPN+i + Pn-1 - 2P, 



N ■ 



(15) 



Formally we can let the number of particles N run from 
—oo to +0O. The matrix s — W can easily be diag- 
onalized. The eigenvalues are s + 2A(1 — cos/c) with 
— TT < A: < TT, and the corresponding eigenvector has com- 
ponents ex-p[ikN]. We can use this to calculate elements 
of the inverse of the matrix s — W. 



1 r „ cos(|iV-M|fc) 

— / dk ; 

27r s + 2A(1 - cos A:) 



1 

ys 



y 



1 



-| \N-M\ 



with 



(16) 



(17) 



This then yields for the probability that a particle ad- 
sorbs and desorbs at a later time without being converted 
the expression 



N,N+1 



(W'rx) 



1 



1 



(18) 



(Here and in the rest of the paper we set s = Wj-x in the 
definition of y.) For the reactivity i?prod we have 

5prod = A[l - fN,N + l{Wr^)] = ^W,^{y - 1). (19) 

We see that the reactivity no longer depends on the load- 
ing 9 if the system becomes large, but only on the reaction 
rate and on A. For small rates of adsorption and desorp- 
tion we again find i?prod ^ A. If the rates of adsorption 
and desorption are large, then 



B 



prod 



(AW,,) 



1/2 



(20) 



for A/Wrx — > oo. The approach to the limit S' — > oo 
becomes very slow when A is large as can be seen in 

Fig. g. 

The procedure for the first-passage problem above can 
also be used to derive the A and B profiles: i.e., the 
distribution of the A's and B's in the chain. We will first 
deal with the question what the probability is that the 
nth particle, counting from site 1, is an A. Then we will 
answer the question that site m is occupied by an A. The 
answers to the same questions about B's follow trivially 
from the ones for the A's. 



We are only interested in the steady state. In this case 
we have detailed balance: i.e., the number of transitions 
per unit time from A'' to + 1 particles is equal to those 
of A^ -I- 1 to N. This means that for each sequence in the 
number of particles 



Nt ^ No 



(21) 



where after a time lapse Ati the number of particles 
changes from A^^ to A^i+i, there is another sequence 



iVo = A^T Nt^ 



AtT-2 



Ati 



A^o. (22) 



The second sequence is the time reversed of the first one. 
Moreover, both sequences are equally likely because of 
the detailed balance. Consequently, the probability dis- 
tribution that the number of particles in the system is 
K at time t = and M with AI < K at time — i for 
the last time equals fMK{t)- This means that if there 
are A^ particles in the system the probability distribu- 
tion that particle n is in the system for a time t equals 
fN-n,N{t)- The probability that that particle has not 
been converted is then fN-n,N{Wrx) following the rea- 
soning after Eq. (||). This probability can be calculated 
using Eq. (|). 

For the profile we need the probability that site m is 
occupied by particle n. This probability is given by 



Pocc(n, A; m, S) = 



N -n 



I 



(23) 



The probability (A^) that site m is occupied by a particle 
that has not been converted is then given by 



(A„ 



in 

E 

n=l 



(24) 



m 

n - 



S — m 

N -n 



lN~n,N 



(M/,x) 



where the first summation averages over the number of 
particles in the system and P^"'' is given by Eq. ( p^ ) 

The expression above can be simplified and interpreted 
more readily for an infinite system. With Eq. (^6|) we 
have 



■,N — 



y+l 



(25) 



for the probability that particle n has not been converted. 
We see that this probability decreases exponentially. 

The probability Pocc{n, N; m, S) becomes Pocc{n, m; 9) 
with 9 = N/S for 5*^00. This limit of the combinato- 
rial factors yields 



Poccin,m] 9) = 



m 

n - 



1 



'(1- 



(26) 
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FIG. 2: The probabilities (An) that a site is occupied by a 
non-converted particle as a function of site index n for rate 
of adsorption and desorption A = 16, loading 6 = 0.2, and 
various system sizes S. 



Substituting this expression and (^5|) in (|2J) yields 

6{x - s)(x 4 



(A, 



26*5)^ 



{x + sY 



(27) 



Again we find an exponential decrease but now with a 
characteristic length of A = l/[ln(y + 1) — ln(y + 1 — 
26*)]. For high loadings {9 1) we find A « y/2. For 
low loading {9 0) we get A « (y + l)/(26'). So the 
characteristic length is larger for low loadings and higher 
rate of adsorption and desorption. 

Figure |^ shows some typical profiles. The straight line 
corresponds to the exponential decrease of Eq. ( p7| ) . The 
result for finite system sizes can be understood from the 
fact that smaller systems are less reactive because there 
are fewer sites at which conversion can take place. As 
the reactivity must be equal to the number of desorbing 
converted particles, Wdcs(Bi) must be smaller for smaller 
S, and because (Bi) = 9 — (Ai) the curves in Fig. || 
must start out at larger values of (Ai). The reactivity 

also equals WrxX)n=i(^")- '^^i^ means that 

must be smaller for smaller S. Because (Ai) is larger for 

smaller S, (A„) must decrease faster. 

Figure ^ confirms the slow convergence of the system 
to the limit S oo. At S' = 100 the values of (A„) drop 
down to 2 • lO^'^. At such low value one would not expect 
an influence of these sites on the kinetics, but one clearly 
sees differences between S = 100 and S" — > oo for all sites. 

We would finally like to comment on three extensions 
of our model: adsorption and desorption at both ends, 
finite diffusion, and other reactions. A system with 2S 



sites and adsorption and desorption at both ends has a 
reactivity that is twice the one of a system with S sites 
and adsorption and desorption at just one end provided 
that S ^ A. DMC simulation show that if S is near or 
smaller than A, then the system open at both ends is 
relatively more reactive. 

Finite diffusion is expected to be equivalent to infinite 
diffusion as long as it is much faster than the other re- 
actions. If this is not the case, then it may become a 
rate determining step. The system then becomes trans- 
port limited. As a consequence the reactivity can only 
go down. We have seen that in our simulations. pO[ 

We expect that our results will change only little when 
we change the details of the reaction or when interaction 
between the particles are included. Our results are a 
consequence of the non-passing of the particles and of the 
fact that particles that stay longer in the system have a 
higher probability of being converted. We therefore also 
do not expect our results to change when we change our 
model of discrete sites to some continuous model. 
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